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Non-perturbative QCD effective charges 
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Using gluon and ghost propagators obtained from Schwinger-Dyson equations (SDEs), we construct the non- 
perturbative effective charge of QCD. We employ two different definitions, which, despite their distinct field- 
theoretic origin, give rise to qualitative comparable results, by virtue of a crucial non-perturbative identity. 
Most importantly, the QCD charge obtained with either definition freezes in the deep infrared, in agreement with 
theoretical and phenomenological expectations. The various theoretical ingredients necessary for this construction 
are reviewed in detail, and some conceptual subtleties are briefly discussed. 
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1. Introduction 

One of the challenges of the QCD is the self- 
consistent and physically meaningful definition 
of an effective charge. This quantity provides 
a continuous interpolation between two physi- 
cally distinct regimes: the deep ultraviolet (UV), 
where perturbation theory works well, and the 
deep infrared (1R), where non-perturbative tech- 
niques, such as lattice or SDEs, must be em- 
ployed. The effective charge depends strongly on 
the detailed dynamics of some of the most fun- 
damental Green's functions of QCD, such as the 
gluon and ghost propagators pQ. 

In this talk we will focus on two characteris- 
tic definitions of the effective charge, frequently 
employed in the literature. Specifically we will 
consider (i) the effective charge of the pinch tech- 
nique (PT) [2 1 3] and (ii) the one obtained from 
the ghost-gluon vertex [3]. 

2. Two non-perturbative effective charges 

We first introduce the notation and the basic 
quantities entering into our study. In the covari- 
ant renormalizable (R$) gauges, the gluon prop- 
agator A MI/ (g) has the form 



A^O?) = -i 



where £ denotes the gauge-fixing parame- 
ter, P^v{q) = g^v — q_iiq_v /q 2 is the usual trans- 
verse projector, and A^ 1 (q 2 ) = q 2 + iH(q 2 ), with 
U MV (q) = P )lv (q)T\(q 2 ) the gluon self-energy. In 
addition, the full ghost propagator D(q 2 ) and 
its dressing function F(q 2 ) are related by 
iF(q 2 ) = q 2 D(q 2 ). The all-order ghost vertex will 
be denoted by F^k, q), with k representing the 
momentum of the gluon and q the one of the anti- 
ghost; at tree-level F^(k,q) = — gv 

An important ingredient for what follows is the 
two-point function A /j;y (g) defined by |H5j 

iA^(?) = A f Hf}D{k + q)^{k)H av {k,q), 

Jk 



9^G{q 2 



q^qu 



L(q 2 ), 



(2) 



(1) 



where A = g 2 Cp Vl Ca is the Casimir eigen- 
value of the adjoint representation, and 
J k = /! 2e (27r)~ d Jd d k, with d the dimension of 
space-time. H fil ,(k, q) is represented in Fig. [I] and 

its tree-level, Hj£) = ig^v In addition, H^(k,q) 
is related to F^k, q) by 

q"H flv (k,q) = -iF tt (k,q). (3) 

2.1. The pinch technique effective charge 

The PT definition of the effective charge relies 
on the construction of an universal (i.e., process- 
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H m (k,q) = Hi a J + 

'/ '••• •o- y 

k + q 

Figure 1. Diagrammatic representation of H. 

independent) effective gluon propagator, which 
captures the running of the QCD (3 function, ex- 
actly as happens with the vacuum polarization in 
the case of QED 31] (See Fig. E2). One impor- 
tant point, explained in detail in the literature, 
is the (all-order) correspondence between the PT 
and the Feynman gauge of the BFM [2 7J. In fact, 
one can generalize the PT construction [2] in such 
a way as to reach diagrammatically any value of 
the gauge fixing parameter of the BFM, and in 
particular the Landau gauge. In what follows we 
will implicitly assume the aforementioned gener- 
alization of the PT, given that the main identity 
we will use to relate the two effective charges is 
valid only in the Landau gauge. 

To fix the ideas, the PT one-loop gluon self- 
energy reads 



A-V) = , 2 



1 



(4) 



where b = 11CU/487T 2 is the first coefficient of the 
QCD /3-function. Due to the Abelian Wis sat- 
isfied by the PT effective Green's functions, the 
new propagator-like quantity A^ 1 (q 2 ) absorbs all 
the RG-logs, exactly as happens in QED with 
the photon self-energy. Then, the renormaliza- 
tion constants of the gauge-coupling and of the 
PT gluon self-energy, defined as 



50 2 ) = Z g 1 (p 2 )g , 

= Z A V)A (<Z 2 ), 



A(g 2 ,/x 2 



(5) 



where the "0" subscript indicates bare quantities, 
satisfy the QED-like relation Z„ 



?-l/2 



There- 



fore, the product 

d (q 2 ) = g 2 A (q 2 ) = 5 V)A(g 2 ,M 2 ) = d(q 2 ), (6) 

forms a RG-invariant (/i-independent) quan- 
tity [2]. For asymptotically large momenta one 
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Figure 2. The universal PT coupling. 




may extract from d(q 2 ) a dimensionless quantity 
by writing, 



d{q 2 



(7) 



where g 2 (q 2 ) is the RG-invariant effective charge 
of QCD; at one-loop 



l + 6. 9 2 ln(?7M 2 ) b\n(q 2 /A 2 QCD ) 



(8) 



where A QCD denotes an RG-invariant mass scale 
of a few hundred MeV. 

Eq. ^ is a non-perturbative relation; there- 
fore it can serve unaltered as the starting point 
for extracting a non-perturbative effective charge, 
provided that one has information on the IR be- 
havior of the PT-BFM gluon propagator A(g 2 ). 
Interestingly enough, non-perturbative informa- 
tion on the conventional gluon propagator A(q 2 ) 
may also be used, by virtue of a general rela- 
tion connecting A(q 2 ) and A(q 2 ). Specifically, a 
formal all-order relation known as "background- 
quantum" identity [8] states that 



A(q 2 



[l + G(q 2 )] 2 A(q 2 ). 



(9) 



Note that, due to its BRST origin, the above re- 
lation must be preserved after renormalization. 
Specifically, denoting by Z A the renormalization 
constant relating the bare and renormalized func- 
tions, Aft" and , through 

«T +A""W) = Z A (p 2 )[g^+A^(q)}, (10) 



then from 
tional relation 

z- 1 = z\ /2 z A , 



and Z g = Z A 



1/2 



follows the addi- 



(11) 



which is useful for the comparison with the cou- 
pling discussed in the following subsection. 

It is now easy to verify, at lowest order, that 
the 1 + G(q 2 ) obtained from Eq. ^ restores the 
j3 function coefficient in front of UV logarithm. 
In that limit [5] 



l + G^ 2 ) 
A-V) 



1 



9C A g 2 



In 



4 48tt 2 \fi 
, , 13C Aff 2 



2 48tt 2 
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(12) 
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Using Eq. © we therefore recover the A _1 (<7 2 ) of 
Eq. (T4|), as we should. 

Then, non-perturbatively, one substitutes into 
Eq. © the 1 + G(q 2 ) and A(q 2 ) obtained from 
either the lattice or SD analysis, to obtain A(q 2 ). 
This latter quantity is the non-perturbative gen- 
eralization of Eq. 0; for the same reasons ex- 
plained above, the combination 



d(q 



g 2 A(q 2 



[1 + G(g 2 )] 2 ' 
is an RG-invariant quantity. 



(13) 



2.2. The ghost-gluon vertex charge 

In principle, a definition for the QCD effective 
charge can be obtained starting from the vari- 
ous QCD vertices; however, in general, such a 
construction involves more than one momentum 
scales, and further assumptions about their values 
need be introduced, in order to express the charge 
as a function of a single variable. The ghost-gluon 
vertex has been particularly popular in this con- 
text, especially in conjunction with Taylor's non- 
renormalization theorem and the corresponding 
kinematics [4]. 

We next define the following renormalization 
constants 



A(q 2 ) = Z^Aoiq 2 ), 
T»(k,q)=Z 1 T»(k,q), 



F(q 2 ) = Z-'F (q 2 ), 
(14) 



Notice that a priori 



Zg' — Z\Z 



-1/2 , 



z g , 



defined as 



does not have to coincide 



with the Z g introduced in ([5]); however, as we 
will see shortly, they do coincide by virtue of the 
basic identity we will derive in next section. 

It turns out that for the so-called Taylor kine- 
matics (vanishing incoming ghost momentum, 
fc M — > —qp), one may impose the additional con- 
dition 



Z 1 = Z g ,zY Z c = l Z-} = Z l J l 2 Z c . (15) 

Thus, the combination 



?(q 2 ) = g /Z A(q 2 )F 2 (q 2 ), 



(16) 



is a RG-invariant (^(-independent) quantity. 
Therefore, for asymptotically large q 2 , in anal- 
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Figure 3. The SDE for the ghost propagator. 



ogy to Eq. one can define an alternative QCD 
running coupling as 



st 2^ ggh(g 2 ) 



(17) 



It is easy to verify that g 2 h (q 2 ) and g 2 (q 2 ) displays 
the same one-loop behavior, since, perturbatively, 
the function 1 + G(q 2 ) is the inverse of the ghost 
dressing function F(q 2 ) (this is due to the general 
identity of Eq.EO]). 

3. An important identity 

In this section, we discuss a non-trivial identity, 
valid only in the Landau gauge, relating the F(q 2 ) 
with the G(q 2 ) and L(q 2 ) of ©. 

The derivation proceeds as follows. First, con- 
sider the standard SD equation for the ghost 
propagator, represented in Fig. [3J and written as 

iD- 1 {q 2 ) = q 2 +i\ f r»A^(k)T»(k,q)D(p),(18) 



where p — k + q. Then, contract both sides of the 
defining equation @ by the combination q^q v to 
get 



{G(q 2 )+L(q 2 )}q 2 =\ q p A<>° (k) q" H„„(k, q)D(p).(l9) 
Jk 

Using the Eq. ((3|) and the transversality of the 
full gluon propagator, we can see that the rhs of 
Eq. (|19p is precisely the integral appearing in the 
ghost SDE (jT8j) . Therefore, in terms of the ghost 
dressing function F(q 2 ), 



l + G(q 2 )+L(q 2 )=F- 1 ( q 2 ). 



(20) 



Eq. (j20)) . derived here from the SDEs, has been 
first obtained in [9 , as a direct consequence of 
the BRST symmetry. 

Let us study the functions G(q 2 ) and L(q 2 ) 
more closely. From Eq. @ we have that (in d 



4 



dimensions) 
G(q 2 ) = 

Hq 2 ) = 



1 



(d-l)ga 
1 



(q 2 A£ - q^A^) , 
2 (dq^A^ - q 2 A») . (21) 



(d- l)q 2 v " * * " 

In order to study the relevant equations fur- 
ther, we will approximate the two vertices, H^ v 
and T^, by their tree- level values, Then, setting 
f(k,q) = (k ■ q) 2 /k 2 q 2 , one may show that pQ 



F-\q 2 ) = l 



+ A /[l- 

k 



f(k,q)]A(k)D(k + q), 



(d - l)G(q 2 ) = A [(d - 2)+/(fc, q)}A(k)D(q + k), 
Jk 



(d - l)L{q 2 ) = A J[l-d f(k, q)]A(k)D(q + k) ,(22) 

It turns out that if F and A are both IR finite, 
Eq. (H2J yields the important result L(0) = pQ. 

Of course, all quantities appearing in Eq. (|22"|) 
are unrenormalized. It is easy to recognize, for 
example, by substituting in the corresponding in- 
tegrals tree-level expressions, that F~ 1 (q 2 ) and 
G(q 2 ) have exactly the same logarithmic depen- 
dence on the UV cutoff, while L(q 2 ) is finite at 
leading order. 

Since the origin of (|20|) is the BRST symme- 
try, it should not be deformed after renormaliza- 
tion. To that end, using the definition of (TIT))) , in 
order to preserve the relation (|20p we must im- 
pose that Z\ = Z c . In addition, by virtue of ([3]), 
and for the same reason, we have that, in the 
Landau gauge, F„(fc, q) and H av {k,q) must be 
renormalized by the same renormalization con- 
stant, namely Z\ [viz. Eq. I|14p]: for the Taylor 
kinematics, we have that Z% = 1 [see Eq. (|15l) ] 
(for some additional subtleties see p]). 

Returning to the effective charges, first of all, 
comparing Eq. ((6]) and Eq. (fT6|) . it is clear that 
g(lj) — g'(n), by virtue of Z\ = Z c . Therefore, 
using Eq. ([9]), one can get a relation between 
the two RG-invariant quantities, r(q 2 ) and d(q 2 ), 
namely 

?{q 2 ) = [\ + G{q 2 )] 2 F 2 {q 2 )d{q 2 ). (23) 

From this last equality follows that a PT (q 2 ) and 
a gh(q 2 ) are related by 
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Figure 4. Numerical solutions for the gluon prop- 
agator. 



After using Eq. (|20|) . we have that 



a PT (<7 2 ) = a gh (<? 2 ) 



1 



L(q 2 



1 + G(<? 2 ) 



(25) 



a gh (q 2 ) = [l + G(g 2 )] 2 F 2 ( g > PT (q 2 ). 



(24) 



Evidently, the two couplings can only coincide at 
two points: (i) at q 2 — 0, where, due to the fact 
that L(0) = 0, we have that a g h(0) = a PT (0), 
and (ii) in the deep UV, where L(q 2 ) approaches 
a constant. 

4. The nonperturbative analysis 

We next turn to the dynamical information re- 
quired for the various ingredients entering into 
the effective charges defined above. To that end, 
we solve numerically the system of SDEs for 
A(q 2 ) F(q 2 ), G{q 2 ) and L{q 2 ) obtained in [5] 

In Figs. J3J and (|5|) we show the results 
for A(q 2 ) and F(q 2 ) renormalized at three 
different points, fi = {4.3, 10, 22} GeV with 
a{n 2 ) = {0.21,0.16,0.13} respectively On the 
right panel we plot the corresponding F(q 2 ) 
renormalized at the same points. Notice that the 
solutions obtained are in qualitative agreement 
with recent results from large- volume lattices [TU] 
where the both quantities, A(q 2 ) and Fq 2 ), reach 
finite (non-vanishing) values in the deep IR. 

The results for 1 + G(q 2 ) and L(q 2 ), renor- 
malized at the same points, are presented in 
Fig. [SI As we can see, the function 1 + G(q 2 ) 
is also IR finite exhibiting a plateau for values 
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Figure 5. Numerical solutions for the ghost dress- 
ing function. 

of q 2 < O.lGeV 2 . In the UV region, we instead 



recover the perturbative behavior (fT2"|) . On the 
other hand, L(q 2 ), Fig. [71 shows a maximum in 
the intermediate momentum region, while, as ex- 
pected, L(0) = 0. 

With all ingredients defined, the first thing one 
can check is that indeed Eq. (fT5)) gives rise to a 
RG-invariant combination. Using the latter def- 
inition, we can combine the different data sets 
for A(q 2 ) and [1 + G(q 2 )] 2 at different renormal- 
ization points, to arrive at the curves shown in 
Fig. [51 Indeed, we see that all curves, for differ- 
ent values of /i, merge one into the other proving 
that the combination d(q 2 ) is independent of the 
renormalization point chosen. 
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Figure 7. L(q 2 ) determined from Eq. (|22 



RGI product d(q 2 }=g 2 £(q ! ) 
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Figure 8. The d(q 2 ) obtained by combining A(q 2 ) 
and [1 + G(g 2 )] 2 according to ([13]) . 



From the dimensionful d(q 2 ) we must now ex- 
tract a dimensionlcss factor, ~g 2 (q 2 ), correspond- 
ing to the running coupling (effective charge). 
Given that A(g 2 ) is IR finite (no more "scal- 
ing"!), the physically meaningful procedure is 
to factor out from d(q 2 ) a massive propagator 
[q 2 +m 2 (q 2 )}-\ 



2, FV 



(26) 



m 2 (q 2 ) ' 

where for the mass we will assume "power-law" 
running [TT], m 2 (q 2 ) = m,Q/(q 2 + to 2 ,). 

Thus, it follows from Eq. ([26]) . that the effective 
charge a PT — ~g 2 (q 2 ) / 4tt is identified as being 



47ra PT ( g 2 ) = [q 2 +m 2 {q 2 )]d(q 2 ), 



(27) 



q [GeV] 



Figure 6. 1 + G{q 2 ) determined from Eq. 
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Figure 9. a gh (q 2 ) vs ap T (q 2 ), for m = 500 MeV. 

Finally we compare numerically the two effec- 
tive charges, a PT (q 2 ) and a g h(q 2 ) in Fig. [HI First, 
we determine a PT (g 2 ) obtained using (|27|) , then 
we obtain a g h(q 2 ) with help of (|25[) and the re- 
sults for 1 + G(q 2 ) and L(q 2 ), Fig. 1 and Fig. [3 

As we can clearly see, both couplings freeze at 
the same finite value, exhibiting a plateau for val- 
ues of q 2 < 0.02 GeV 2 , while in the UV both show 
the expected perturbative behavior. They differ 
only slightly in the intermediate region where the 
values of L(q 2 ) are appreciable. 



5. Conclusions 

In this talk we have presented a comparison 
between two QCD effective charge, apT(q 2 ) and 
a g h(<z 2 ), obtained from the PT and the ghost- 
gluon vertex, respectively. 

Despite their distinct theoretical origin, due to 
a fundamental identity relating the various ingre- 
dients entering into their definitions, the two ef- 
fective charges are almost identical in the entire 
range of physical momenta. In fact, the coincide 
exactly in the deep infrared, where they freeze at 
a common finite value, signaling the appearance 
of IR fixed point in QCD [12] , also required from 
a variety of phenomenological studies [13] . 
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